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Abstract 

We investigate the composition and shape distribution of silicate dust grains in the interstellar medium. The effects of the amount of magnesium 
and iron in the silicate lattice are studied in detail. We fit the spectral shape of the interstellar 10pm extinction feature as observed towards 
the galactic center using various particle shapes and dust materials. We use very irregularly shaped coated and non-coated porous Gaussian 
Random Field particles as well as a statistical approach to model shape effects. For the dust materials we use amorphous and crystalline silicates 
with various composition as well as silicon carbide (SiC). The results of our analysis of the 10/jm feature are used to compute the shape of 
the 20 fjm silicate feature and to compare this with observations of this feature towards the galactic center. By using realistic particle shapes 
to fit the interstellar extinction spectrum we are, for the first time, able to derive the magnesium fraction in interstellar silicates. We find that 
the interstellar silicates are highly magnesium rich (Mg/(Fe + Mg) > 0.9) and that the stoichiometry lies between pyroxene and olivine type 
silicates (O/Si ~ 3.5). This composition is not consistent with that of the glassy material found in GEMS in interplanetary dust particles 
indicating that the amorphous silicates found in the Solar system are, in general, not unprocessed remnants from the interstellar medium. Also, 
we find that a significant fraction of silicon carbide (~3%) is present in the interstellar dust grains. We discuss the implications of our results 
for the formation and evolutionary history of cometary and circumstellar dust. We argue that the fact that crystalline silicates in cometary and 
circumstellar grains are almost purely magnesium silicates is a natural consequence of our findings that the amorphous silicates from which 
they were formed were already magnesium rich. 
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1. Introduction 

Understanding the composition of silicate dust in the 
Interstellar Medium (ISM) is crucial for studies of circumstel- 
lar and cometary dust and for the processing of dust in in- 
terstellar space. For example, in order to understand the way 
in which dust is modified in protoplanetary disks during the 
process of disk dissipation and planet formation, it is impor- 
tant to know the original composition of the material that en- 
tered the disk. It is reasonable to assume that this material 
had a composition similar to that found in the diffuse inter- 
stellar medium. The amount of iron and magnesium in the sil- 
icates contains valuable information on the origin of the grains 
and is an important probe of the processing history. When 
the silicates are crystalline, i.e. they have a long range or- 
der in the lattice, the magnesium content can be readily ob- 
tained from the spectral positions of sharp solid state emis- 
sion features. In this way it has been shown that cometary 
and circumstellar crystalline silicates are highly magnesium 
rich, i.e. Mg/(Fe + Mg) > 0.9 (see e.g. IjageretalJ, U998 : 



Fabi an et all I200U iKoike et al.L I2003I) . For amorphous sili- 
cates, i.e. silicates with a disordered lattice structure, the spec- 
tral features are less pronounced which makes it harder to de- 
termine the composition from spectroscopy. In particular, it is 
hard to disentangle effects of grain shape and composition on 
the spectral profile. In previous studies the interstellar silicates 
have been modeled using amorphous silicates containing ap- 
proximately equal amounts of iron and magnesium (see e.g . 
Draine & Leelll984tlLi & Drainell200ltlKemper et alll2004) . 



By inclu ding amorphous sili cates with slightly different com- 
positions Kemper et alj (12004b concluded that the composition 
is indeed most likely Mg/(Fe + Mg) w 0.5 - 0.6. If the grain 
properties derived in these studies are correct it is difficult to 
understand the high magnesium content of crystalline silicates 
in protoplanetary disks and solar system comets. The reason 
is that the majority of crystalline silicates is expected to be 
formed by thermal annealing of amorphous silicates, preserv- 
ing the composition of the amorphous silicates, while only a 
small fraction is expected to originate from gas-phase conden- 
sation, producing magnesium rich crystals dGaill 120041) . 
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Figure 1. Typical examples of porous Gaussian Random Field particles generated using the method described in the text. Every 
small cube seen in this figure represents a volume element according to the discretization of space used for the construction of 
the Gaussian random field. The grains shown here are typical examples from the ensemble of 38 porous GRF particles used to 
compute the extinction spectra. 



In the above mentioned studies the 10 //m ISM extinction 
feature has been modeled assuming that the optical properties 
of the grains can be represented by those of perfect homoge- 
neous spheres. This allows for fast and simple computations. 
The fit to the 10 //m ISM extinction feature using homogeneous 
spheres is surprisingly good, and it is therefore tempting to as- 
sume that the derived dust properties are correct. There are, 
however, several indications that the assumption of homoge- 
neous spheres cannot be adequate. First there is the interstellar 
polarization. A homogeneous sphere of optically inactive ma- 
terial, like silicates, cannot produce linear polarization through 
extinction. Since the lO/vm silicate feature is clearly visible in 
the spectr al dependence of the degree of interstella r linea r po- 
larization dLee & DraineLll985l:lAifken et al.L[l986Lll989h . the 
silicates have to contribute to it. The most natural explanation 
for this is extinction by aligned, nonspherical silicate grains. 
The interstellar polarization alone, however, does not preclude 
using the spherical grain model for explaining the extinction 
spectrum. One might assume that spherical grains can be used 
to model the average extinction properties of slightly elongated 
particles. In addition, it appears that only very modest devia- 
tions from a perfect sphere are needed to produce the required 
polarization (an aspect ratio of ~ 2 seems sufficie nt, see e.g. 
Lee & Draine . 1985 : Henning & S tognienkol 1993h . However, 
it has been shown recently that the 20 //m silicate extinction 
band cannot be accurat ely reproduced by usin g homogeneous 
spherical dust grains dChiar & Tielensl 120061) . In addition it 
is clear from several studies that effects of particle shape on 
the extinction spectrum of small particles can be very large. 
Therefore, a detailed re-analysis of the 10/im ISM extinction 
feature using more realistic grain models is needed. 

In this paper we study the effects of grain non-sphericity 
on the derived composition and lattice structure of interstellar 
grains. We do this by computing the optical properties of an 
ensemble of so-called porous Gaussian Random Field (GRF) 
particles and by applying a statistical approach to account for 
particle shape effects. For this study the exact choice of the 
shapes of the nonspherical particles is not crucial since it is 
shown by iMin et al.1 ( 12003b that the optical properties of var- 
ious classes of nonspherical particles are quite similar, and it 



is mainly the perfect symmetry of a homogeneous sphere that 
causes a different spectral signature. From our analysis of the 
10 /urn spectral extinction feature we can compute the expected 
shape of the 20 fim feature. We compare this to the ob served 
extinction spectrum obtained bv lChiar & Tielensl (120061) . 

The structure of the paper is as follows. In section [2] we 
explain in some detail the effects of particle shape and compo- 
sition on the shape of the 10 pun spectral profile. The procedure 
to fit the 10 fim feature and the results of this fit are presented 
in section [3] Implications of our results are discussed in sec- 
tion @] Finally, in section [5] we summarize the conclusions of 
the paper. 



2. Extinction spectra 

2.1. Effects of grain shape 

The shape of a dust grain has a large influence on the shape of 
the absorption and ext inction spectrum it causes. This effect has 



been described by e.g.lBohren & Huffman ( 1983l):lFabian et al . 
d200ll) ; lMin et al.ld2003h : IVoshchinnikov et al.ld2006l) . In many 



cases one needs to consider irregularly shaped dust grains 
in order to reproduce the infrared spectra observed from cir- 
cumstellar cometary or laboratory particles correctly (see e.g . 



Honv et al., 2002; Bou wman et all l2003t iMoreno et aU I2003 
IMin et all I2003I) . It can be shown that the absorption spec- 
trum caused by homogeneous spherical particles is very differ- 
ent from that caused by other particle shapes. This difference 
is much larger than the differences du e to various nonspheri- 



cal particles shapes (IMin et a l.. 2003). The perfect symmetry 



of a homogeneous sphere causes a resonance with a spectral 
shape and position very different from resonances observed in 
the spectra of all other particle shapes. This emphasizes why it 
is puzzling that the interstellar lO/zm extinction feature can be 
fitted almost perfectly using homogeneous amorphous silicate 
spheres. In this paper we will consider three types of particle 
shapes, homogeneous spheres, porous Gaussian Random Field 
(GRF) particles and a statistical ensemble of simple particle 
shapes to represent irregularly shaped particles. 
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Figure 2. A porous GRF particle without (left) and with (right) 
a carbon mantle. Both particles are created from the same 
Gaussian Random Field. The dark grey cubes in the right pic- 
ture represent the amorphous carbon mantle which takes up 
30% of the material of the grain. 



2.1.1. Gaussian Random Field particles 

Several attempts have been made to properly model irregularly 
shaped particles. One of the most widely used classes of ir- 
regularly shaped particles is perh aps that of Gaussian Random 
Spheres (Mui nonen et al. . 1996). These particle shapes have 
the advantage that they are easily generated numerically. 



However, the resulting particles do not appear to resemble ac- 
curately realistic particle shapes thought to be representative 
for those in the ISM li ke, for example, I nterplanetary Dust 



Particles (IDPs; see e.g. Warr en et all 119941) . A more promis 



ing class of particles in this respect mig ht be the Gaussian 
Random Field (GRF) particles (se e e.g. iGrvnko & Shkuratovi 
2003 ; Shkuratov & Grvnkol 2005 ). This class of particles cov- 
ers a wide range of very irregular shapes. In addition, the par- 
ticles are also relatively easily generated numerically. 

In short, a GRF particle is constructed as follows. First, 
space is divided into small volume elements. Every volume el- 
ement is then given a value which is correlated with the values 
of surrounding volume elements according to a lognormal dis- 
tribution. The resulting three dimensional field is called a GRF. 
A threshold level is introduced (in our case 0.5 times the max- 
imum value). All volume elements with a value of the GRF 
above the threshold are inside the particle, while all other vol- 
ume elements are considered vacuum. In this way an ensemble 
of particles is created. Since we considered the resulting par- 
ticles to be still rather smooth, we developed a way to make 
the particles more irregular. This was done by creating a sec- 
ond GRF with a much smaller correlation length. Where this 
second field has a value below a certain threshold, we place 
vacuum voids. In this way, irregularly shaped voids are created 
inside the particles, making them porous, and on the edge of the 
particles, making them more rough. We will refer to these par- 
ticles as porous GRF particles. Some pictures of typical exam- 
ples of porous GRF particles are shown in Fig.Q] Every small 
cube seen in this figure is a volume element according to the 
discretization of space chosen for the construction of the GRF. 
We construct an ensemble of 38 porous GRF particles and av- 
erage the extinction spectra over this ensemble. For details on 




A (ftm) 

Figure 3. The mass extinction coefficient as a function of wave- 
length in the 10 //m region for amorphous silicate particles 
with an olivine type stoichiometry and x = Mg/(Fe + Mg) = 
0.5 (i.e. MgFeSiO/0 with various shape distributions. The 
feature as obtained using the widely employed Continuous 
Distribution of Ellipsoids (CDE; grey line) is shown for com- 
parison. 



the algorithm to construct porous Gaussian Random Field par- 
ticles, see Appendix lAl 

In order to study the effect of mixed carbon/silicate grains 
we also consider coated porous GRF particles. These are cre- 
ated in such a way that a fixed abundance of 30% of the outer 
layer of the particle consists of amorphous carbon. This can be 
done easily using porous GRF particles by taking two thresh- 
olds in the GRF from which the particle is created. For points in 
the grain where the value of the GRF is > 0.5 and smaller than 
the second threshold, we have the carbon mantle. For points 
in the grain where the GRF is larger than the second threshold, 
we have the core of the grain. The value of the second threshold 
is chosen such that the mantle contains 30% of the total grain 
material volume. A picture of a coated porous GRF particle is 
shown in Fig. [2] 

2.1 .2. Distribution of Hollow Spheres 

When modeling the optical properties of irregularly shaped par- 
ticles in astrophysical environments we have to keep in mind 
that we always observe the average optical properties of an 
ensemble of particles with various shapes. Therefore, we are 
not interested in the specific characteristic optical properties 
of a single particle, but rather in the average optical proper- 
ties of the entire ensemble. Once we realize this, it is a logical 
step to use a statistical approach to compute the average op- 
tical properties of an ensemble of irregularly shaped particles. 
In this statistical approach we assume that the optical proper- 
ties of an ensemble of irregularly shaped particles can be mod- 
eled in a statistical way by using the average optical properties 
of a distribution of simple particle shapes the individual parti- 
cles of which need not necessarily be realistically shaped. For 
particles in the Rayleigh domain, i.e. particles much smaller 
than the wavelength both inside and outside the particle, it was 
shown analytically that for all ensembles of particle shapes an 
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ensemble of spheroidal shape s exists with iden tical absorption 
and extinction cross sections dMin et al.ll2006l) . One advantage 
of using the statistical approach is that the computation time 
is limited since it avoids computations of large ensembles of 
complex shaped particles. Another advantage is that by using 
simple particle shapes, the number of parameters used to de- 
scribe the particles is limited. This helps determining the main 
characteristics of the particles and avoids over-interpretation. 

A particularly successful implementation of the statistical 
approach is given by the Distribution of Hollow Spheres (DHS; 



Min et al. , 2003, 2005a). With the DHS it is possible to repro- 



duce measured and observed spectra of small particles while 
computations can be performed easily and fast for ar bitrary par- 



ticle si zes (for an application of this method see e.g. lMin et al. 



2005b). The DHS considers hollow spherical particles, aver- 
aging the optical properties uniformly over the fraction of the 
total volume occupied by the central vacuum inclusion, /, over 
the range < / < / max while keeping the material volume 
of the particles constant. The value of / max reflects the degree 
of irregularity of the particles considered. We would again like 
to stress here that by taking this approach we do not consider 
the real particles to be hollow spherical shells. Rather, it is as- 
sumed that the optical properties of an ensemble of realisti- 
cally shaped, irregular particles can be represented by the aver- 
age properties of an ensemble of hollow spheres in a statistical 



2.1 .3. The amorphous silicate feature 

The interstellar silicates causing the lO^um silicate extinc- 
tion feature are mostly amorphous. iKemper et al. I d2004l2005b 
placed an upper limit on the amount of crystalline silicates of 
2.2%. In addition, the particles in the interstellar medium are 
very small (< 0.3 //m). At A — 10/im this implies that the par- 
ticles are much smaller than the wavelength both inside and 
outside the particle, i.e. they are in the Rayleigh domain. For 
particles in the Rayleigh domain the actual particle size is not 
important for the shape of the absorption and extinction spec- 
tra. In order to compute the spectra of por ous GRF particles 
we employed the method of iMin et al. (2006) which allows for 
rapid computation of absorption and extinction spectra of arbi- 
trarily shaped particles in the Ra yleigh dom ain. For the spher- 
ical particles we use Mie theory (Mie, I1908I) . To compute the 
optical properties of the distributi on of hollow spheres we u se 
a simple extension of Mie theory ( Toon & Ackermanl 198 lh . 

Fig. [3] shows the mass extinction coefficient, k, as a func- 
tion of wavelength for amorphous silicate grains with an olivine 
type stoichiometry and an equal amount of magnesium and iron 
(see also section |2~2| | for homogeneous spherical particles, for 
an ensemble of porous GRF particles, for two distributions of 
hollow spheres with / max = 0.5 and 1.0, and for a Continuous 
Distribution of Ellipsoids (CDE). It is clear that the shape and 
position of the feature has a strong dependence on grain shape. 
In general the spectral extinction features caused by irregularly 
shaped particles are much broader and shifted towards the red 
with respect t o those cau sed by homogeneous spherical parti- 
cles (see also IMin et all 120031) . This effect becomes stronger 
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Figure 4. The mass extinction coefficient as a function of wave- 
length in the 10 /im region for porous GRF particles of various 
compositions. 

when deviating more from the perfect homogeneous sphere, as 
can be seen when increasing / max . 

2.2. Effects of grain composition 

Although silicates with various compositions all display a spec- 
tral feature in the 10 /mi region due to the Si-O stretching 
mode, there are differences in the spectral appearance that al- 
low to constrain the composition of the silicates. In our anal- 
ysis we consider amorphous silicates with various composi- 
tions, the crystalline silicates forsterite and enstatite, and sil- 
icon carbide (SiC). Although we expect that carbonaceous 
species and metallic iron are important constituents of inter- 
stellar dust, we have no constraints on their abundances from 
spectroscopy around 10 fim since in this wavelength range they 
have a smooth, featureless extinction spectrum. Therefore, they 
are not detectable in our analysis. The various components we 
use will be discussed below. 

Silicates are composed of linked Si04 tetrahedra with 
cations inbetween (e.g. Mg, Fe, Na, Al). A silicate can 
have various stoichiometries depending on the number of 
oxygen atoms shared between the Si04 tetrahedra. Silicates 
with an olivine type stoichiometry share no oxygen atoms 
and the bonds are formed through the cations. For mag- 
nesium/iron olivines this results in a chemical composition 
Mg2xFe2-2jcSi04, where < x < 1 determines the magne- 
sium over iron ratio, x = Mg/(Fe + Mg). Silicates with a py- 
roxene type stoichiometry share one oxygen atom, resulting in 
linked chains of tetrahedra, which are bound by cations, result- 
ing in a composition Mg v Fei_ x Si03. When all oxygen atoms 
are shared between the tetrahedra we get silica, Si02- Besides 
magnesium/iron silicates we also consider a sodium/aluminum 
silicate with a pyroxene stoichiometry and equal amounts of 
Na and Al, i.e. Nao.sAlo.sSiOs. Since amorphous silicates have 
a disordered lattice structure, we often find silicates in na- 
ture with a composition between the above mentioned types. 
Mixing the silicates mentioned above results in an average 
composition Mg tl Fe X2 Na V3 Al x ,SiO A - 4 , where X4 — x\ + x^ + 
2xj, + 2. We will use X4 = O/Si as a measure for the average 
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Name 


Composition 


Lattice Structure 


Ref. 


Olivine (x = 0.5) 


MgFeSi0 4 


Amorphous 


[1] 


Pyroxene (x = 0.5) 


MgFeSi 2 6 


Amorphous 


[1] 


Olivine (x = 1) 


Mg 2 Si0 4 


Amorphous 


[2] 


Pyroxene (x = 1) 


MgSi0 3 


Amorphous 


[1J 


Na/Al Pyroxene 


NaAlSi 2 6 


Amorphous 


[3] 


Silica 


Si0 2 


Amorphous 


[4] 


Forsterite 


Mg 2 Si0 4 


Crystalline 


[5] 


Enstatite 


MgSi0 3 


Crystalline 


[6] 


Silicon Carbide 


SiC 


Crystalline 


[7] 



Table 1. The materials used in the fitting procedure. In the 
case of the amorphous silicates we use the names 'olivine' or 
'pyroxene' to indicate the average silicate stoichiometry. For 
the coating of the coat ed GRF partic l es we use amorphous 
carbon as measured by Preibisch et al. (1993). The references 
in the last column refer to: [1] Dorschner et al. I (fl995l) . [2] 

1998h. 



Pfenning & Stognienkol d 19961) [3] 



,__ _ , , Mutschke et al 

Soitzer & Kleinmanl dl960h. 151 IServoin & Pirio 
Jager et all d 19981) . and m lLaor & Draind(ll993l 



P dl998h 
(1973). 



14] 
16] 



stoichiometry. For olivines x\ = 4, for pyroxenes x 4 = 3 and 
for silica x\ — 2. 

In general the 10 pun feature for amorphous silicate with an 
olivine type stoichiometry is broader and peaks at longer wave- 
lengths than that of amorphous silicates with a pyroxene type 
stoichiometry. For magnesium/iron silicates with much magne- 
sium, i.e. high values of x, the lOyum feature is shifted to shorter 
wavelengths compared to iron rich silicates, i.e. small values of 
x. The effect of composition on the shape of the lOfim silicate 
feature is illustrated in Fig.|4]for porous GRF particles. 

In our analysis we also include crystalline silicates. As we 
mentioned earlier, the crystalline silicates found in various as- 
tronomical environments are all magnesium rich, i.e. x « 1 . We 
include the two most commonly found crystals namely crys- 
talline forterite (Mg2Si0 4 ) and crystalline enstatite (MgSiC>3) 
in our analysis. Small crystalline silicate grains display very 
strong spectral resonances that are easily detectable in emis- 
sion or extinction spectra. 

Silicon carbide (SiC) is a dust material pr oduced in car- 
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Figure 5. The normalized extinction profile in various lines 
of sight. A straight continuum is subtracted such that at 8.3 
and 12.3 p/m the extinction is 0. The features are normalized 
such that the maximum extinction equals unity. In the upper 
panel we overlay the extinction features towards Sgr A*, GCS3, 
GCS4 and WR98a. In the three panels below this we show the 
latter three extinction profiles minus the profile observed to- 
wards Sgr A*. 



dust grains can, in principle, be determined. In practice, the 
spectrum of the source is often not known, and it is only possi- 
ble to derive the optical depth in the silicate feat ure itself. This 
was do ne for the line of sight towards Sgr A* by lKemper et alJ 
and we will use their results in our analysis. Although 
we restrict ourselves to this line of sight, the extinction fea- 
tures in other lines of sight towards the galactic center look 
very similar (see Fig. |5J cf. Kemper et al.1 120041) . We use the 
feature towards Sgr A* because it has the best signal to noise. 



bon rich AGB stars (see e.g. Speck et al., 1997). Furthermore, 



presolar SiC grains have been fo und in meteorites and IDPs 3.2. Fitting procedure and error analysis 



(see e.g. lBernatowicz et all 119871) . SiC displays a strong spec- 



tral resonance around 1 1 fim. Therefore, we include the possi- 
bility of SiC in our analysis. 

An overview of all dust materials included in our analysis 
with references to the laboratory data used to compute the ex- 
tinction spectra can be found in Table Q] 

3. Interstellar silicate grains 

3.1. Observations 

The optical depth as a function of wavelength for interstellar 
dust grains can be determined by observing a source located be- 
hind significant amounts of interstellar dust, for example near 
the galactic center. When the spectrum of the observed source 
is known the optical depth at each wavelength point, t a , of the 



We fit the interstellar 10 /-im extinction feature by making a lin- 
ear least square fit to the optical depth in the feature using the 
extinction spectra computed for the various dust components 
listed in Table [TJ In order to correct for the continuum opacity 
in the computed extinction spectra of the dust species we sub- 
tract a straight line from the fit. We do this by introducing two 
additional fit parameters to the fitting procedure, the slope and 
offset of this line. In this way we treat the observations and the 
opacities in a consistent manner. We would like to stress here 
that by this procedure we fit the optical depth in the feature, and 
not the total optical depth caused by the dust species. In order 
to get an idea of the errors we make by assuming a straight line 
continuum, we have also made fits using different power law 
continua. This only had a very small influence on the results 
and did not affect any of our conclusions. 
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To avoid fits with an unlikely dust composition due to un- 
realistic combinations of elemental abundances we fix the ele- 
mental abund^nceof Mgrelative to Si. Similar constraints were 
used bv lZubko et al.l ( 120041) to fit the interstellar extinction pro- 
file. To do this we need a good estimate for the abundances 
of Mg and Si in the ISM. Zubko et al.l (120041) use the average 
abundances of F-, G- or B-type stars to constrain their model. 
We use the Mg/Si ratio as derived by Ueda et all ( 20051) from 
X-ray spectroscopy. Using their technique they probe simulta- 
neously the elements in the gas phase and the solid phase of the 
ISM. They derive Mg/Si = 1.22 + 0.14, where the error reflects 
a 90% confidence limit (which corresponds to ~ 1.65cr). We 
use a lcr error implying Mg/Si = 1.22 + 0.09. We assume all 
Mg and Si t o be in the solid phase consistent with the results of 



Ueda et all (2005). 



In order to get an idea of the error on the parameters fit- 
ted to the interstellar 10 /mi extinction feature we generated 
1000 different synthetic observations. For the extinction spec- 
trum we did this by adding noise to the spectrum according 
to the error on the spectrum at each w avelength poin t , i.e. (T a- 
The error on the spectrum as derived bv lKemper et al.l d2004l) is 
ve ry small. The large va lue of x 2 obtained from the spectral fit 
bv lKemper et al.l ((2004) already indicates this. This is probably 
due to deficiencies in the model. To account for this underesti- 
mate of the error we scaled the o~a by a factor y/x* and used this 
scaled error to construct the 1000 synthetic spectra to be fitted. 
In this way when the fit to the spectrum is poor, i.e. a high value 
of x 1 is obtained, the errors on the derived fit parameters will 
also be large. This reflects that a high value ofx 2 indicates that 
the model does not accurately describe the observations. 

For each spectrum to be fitted we also constructed a syn- 
thetic abundance constraint using Mg/Si = 1.22 + 0.09. So we 
employ a slightly different abundance constraint for each syn- 
thetic spectrum. For each of the 1000 combinations of spectra 
and abundance constraints we applied our fitting procedure to 
obtain the abundances of the various dust species. This results 
in a distribution of abundances for all components. This distri- 
bution can be used to compute the average fit parameters along 
with an estimate of the standard deviations on these parameters. 

The errors obtained by the above described procedure are 
the errors on the abundances assuming the dust component is 
present. A better measure for the significance of the detection 
of a dust component can be obtained by performing a so-called 
F-test. This test determines the significance of a component in 
the fit by measuring the decrease of the^ 2 when the component 
is added. A simple F-test, as usually employed, is only valid un- 
der strict assumptions of the distribution of x 2 under variation 
of the observations. However, fitting the 1000 synth etic spec- 
tra al lows us to create the F-distribution exactly (see Che rnickL 
19991) . For each component we also make a fit to all the spec- 



tra leaving out this component. This gives two distributions of 
X 2 , one with and one without the specific component. The dis- 
tribution of the ratios of elements from these two distributions 
is called the F-distribution. By determining the mean and the 
standard deviation, cr, of this distribution, we can determine 
how far (in terms of a multiple of cr) the mean lies from unity. 
Since a mean value of unity implies that the fit is insensitive for 
this specific component, the deviation from unity gives the sig- 



6. 



nificance of the detection. For more details we refer to Chernickl 
(1 19991) . In order to minimize the number of parameters in the 
fit and in this way avoid over interpretation of the observations, 
we iteratively remove the least significant component until all 
components are detected with at least lcr. 
To summarize the fitting procedure: 

1 . Make a linear least squares fit to the observed Tf eature using 
the abundance const raint and the errors o n the spectrum, 
cr A , as determined bv lKemper et al.l (12.004ft . This yields the 
reported^ 2 . 

2. Construct cr^ so that the best fit using these errors would 

have a^- 2 of unity, i.e. a J l = cr A tJx 2 - 

3. Use these cr' A and the errors on the abundance constraint to 
construct 1000 combinations of synthetic spectra and abun- 
dances by adding Gaussian noise to the observations. 

4. Make linear least squares fits to all these synthetic data sets. 

5. For all components in the fitting procedure make linear 
least squares fits to all synthetic data sets leaving out this 
component. 

Use the fit parameters of all 1000 fits to obtain the average 
abundances and use the x 2 distributions obtained in steps 4 
& 5 to construct the F-distributions. The F-distributions are 
used to derive the significances of all dust components. 
If there are dust components detected with a significance 
less than lcr, remove the least significant component and 
restart the fitting procedure from point 1. Iterate until all 
components left in the procedure are detected with at least 

lcr. 



3.3. Results from the 10 micron feature 

The best fits to the 10 pm feature using homogeneous spheres, 
porous GRF particles, coated porous GRF particles, and a 
Distribution of Hollow Spheres are shown in Fig. [6] For the 
DHS we determined that the best fit was obtained by taking 
/max = 0.7, similar to the best fit valu e found in the analysis of 
comet Hale-Bopp dMin et al.[l2005bl) . The derived composition 
for each of these shape distributions is listed in Table[2] The^f 2 
of these fits show that the feature is better represented by using 
irregularly shaped particles than by using homogeneous spher- 
ical particles. 

The resulting dust composition using irregularly shaped 
particles is very different from that obtained when using ho- 
mogeneous spheres. Since it seems unlikely that the interstel- 
lar silicates are perfect homogeneous spheres, and there are no 
other particle shapes that produce a spectrum similar to that 
of perfect homogeneous spheres, this model must be rejected. 
In addition, the fit using homogeneous spheres has a reduced 
X 2 which is almost six times as high as those obtained using 
other particle shapes. From the analysis we draw the following 
conclusions. 

i) The interstellar silicate grains contain much more mag- 
nesium than previously assumed. From Table [2] it can be de- 
duced that when using the porous GRF particles or the DHS, 
we derive that the magnesium content Mg/(Fe + Mg) = 0.91 
and 0.93 respectively. Using homogeneous spheres we derive 
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Figure 6. The wavelength depende nce of the extinction of the best fit models (black curves) together with the extinction feature 
as obtained by lKemper et alJ (|2004) (gray curves). Here, Tf eature denotes the optical depth in the lO/vm feature, and r res idue denotes 
the observed optical depth minus the modeled optical depth. 
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DHS (/ max 
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x 2 




65.3 




12.6 




11.5 




10.1 




Silicates 


Olivine (A; x = 0.5) 


MgFeSi0 4 


57.9 


(6.2a-) 


16.9 


(1.9a-) 


9.7 


(l.lo-) 


13.8 


(2.0a-) 


Pyroxene (A; x = 0.5) 


MgFeSi 2 O s 


12.9 


{1.6a) 














Olivine (A; x = 1) 


Mg 2 Si0 4 


28.9 


(1.0a-) 


47.2 


(3.5a-) 


48.2 


(3.5a-) 


38.3 


(3.8a-) 


Pyroxene (A; x = 1) 


MgSiOj 






24.0 


(6.1a-) 


31.4 


(6.8a-) 


42.9 


(6.7a-) 


Na/Al Pyroxene (A) 


NaAlSi 2 6 






6.8 


(3.2a-) 


5.3 


(2.7a-) 


1.8 


(1.1a-) 


Silica (A) 


SiO z 


0.3 


(1.1a-) 














Forsterite (C) 


Mg 2 Si0 4 






1.5 


(4.7a-) 


1.3 


(3.6a-) 


0.6 


(3.2a-) 


Other 


Silicon Carbide (C) 


SiC 






3.6 


(3.7a-) 


4.2 


(4.3a-) 


2.6 


(5.3a-) 


The average silicate composition implied by the above abundances 


M 


Mg/Si 


1.22 




1.37 




1.39 




1.32 




X2 


Fe/Si 


0.60 




0.13 




0.07 




0.10 




X 3 


(Na, Al)/Si 


0.00 




0.05 




0.03 




0.01 




Xn 


O/Si 


3.82 




3.59 




3.53 




3.45 





Table 2. The composition in terms of mass fractions in percent together with the significance of the detection of the dust compo- 
nents for the best fit models using various particle shapes. A '-' denotes that the dust component was detected with a significance 
less than lcr and was thus removed from the fitting procedure. Enstatite was not found to be significant using any of the four 
grain shapes and is therefore omitted from the table. The (A, C) behind the name in the first column denotes if the material is 
amorphous or crystalline. The names 'olivine' and 'pyroxene' are used for the amorphous silicates to refer to the average stoi- 
chiometry. The x\,X2, x%, X4 refer to the average silicate composition, Mg Al Fe A2 Na X3 Al l3 SiO X4 implied by the abundances of all 
silicates. 



Mg/(Fe + Mg) = 0.67. The fit using coated GRF particles even ii) We confirm the findings by iKemper et alJ ([2004, 2005) 
results in Mg/(Fe + Mg) = 0.95. that the fraction of crystalline silicates (i.e. crystalline forsterite 

and enstatite) in the ISM is very low (~ 0.6 - 1.5%). 
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three different sets of refractive indices. First we considered the 



A (/jm) 

Figure 7. The mass extinction coefficient as a function of wave- 
length in the 10 //m region for silicon carbide (SiC) particles 
with various shape distributions. 



iii) We find a significant fraction (~2.6-4.2%) of silicon car- 
bide (SiC). This implie s that ~9-12% of the available Si is in 
these SiC grains. Kemp er et al. (2004) placed an upper limit of 
0.1% on the abundance of homogeneous spherical SiC grains 
in the ISM from their analysis of the 10 pun spectrum. However, 
the position and shape of the SiC resonance in th e 10 /mi region 
is ver y sensitive to the particle shape (see e.g. lAndersen et all 
20061) " and size (see lSpeck et aUl2005l) Whereas homogeneous 
spherical SiC grains show a sharp resonance at 10.6 /mi, ir- 
regularly shaped SiC particles show a much b roader spec 
tral fe ature peaking near 1 1.25 fim (see Fig. [7] or Speck et al 



20051) . The observed 10/mi extinction spectrum (Fig. [6]) has a 
shoulder around 1 1 fim which is fitted correctly using(coated) 
porous GRF or DHS particles but which is not fitted at all us- 
ing homogeneous spheres. This is due to the presence of ir- 
regularly shaped SiC grains. The presence of SiC in the ISM 
is not unexpected since presolar Si C grains have been foun d 



in dust grains in the Solar system (Ber natowicz et all 119871) 



These SiC grains found in IDPs are indeed irregularly shaped. 
Implications of the presence of SiC in the diffuse ISM will be 
discussed in a separate study. 

iv) We find that the average stoichiometry of the silicates 
corresponds to Si/O ~ 3.5 when applying irregularly shaped 
particles, i.e. between an olivine and a pyroxene type stoi- 
chiometry. When using homogeneous spherical particles we 
find a st oichiometry more w eighted towards olivine type sili- 
cates (cf. lKemper et al. ■ l2004h . The implications for these find- 
ings will be discussed in section|4] 

3.3.1 . The SiC lattice structure 

Silicon carbide can ex ist in many different crystal structures 



Mutschke et al. (1999) have studied the dependence of the 



spectral appearance of the 1 1 /im resonance on crystal struc- 
ture. They concluded that when the crystal is pure, the actual 
crystal structure is not very important. However, they did note 
a large difference when the crystal structure is destroyed and 
an amorphous SiC structure is left. To test the influence of the 
optical data taken to model the SiC resonance we considered 



data as shown before dLaor & Drainei 119931) . Second we take 
a rather pure crystal QS -SiC with y - 10 cm "" 1 ) according to 
the equations given bv lMutschkeetalJj l999). Third we con- 
sider amorphous SiC as measured by [Mutschke et al.l (1999). 
The results of the different fits are shown in Table [3] From this 
analysis it seems that the most likely a ppearance of SiC in th e 
ISM is closest to the measurements by Laor & Drainei dl993l) . 
The abundances of the other dust components are not signifi- 
cantly different from those presented in Table [2] 

3.4. Comparison with the 20 micron silicate feature 

Apart from the frequently used 10 /mi feature, silicates also dis- 
play a spectral signature that peaks near 20 fim due to the O-Si- 
O bending mode. It is often referred to as the 18 /mi feature, 
although its spectral position varies with silicate composition 
and particle shape. We will refer to it as the 20 fim feature. The 
spectral shape of this feature is very sensitive to particle shape 
and composition. Because it is rather broad, it is very hard 
to extract the 20 /mi extinction profile accurately from mea- 
sured infrared spectra. Only if the spectrum of the emission 
source behind the column of extincting dust is known rather 
accurately, a reliable extracti on of the extinction profi le can be 
made. This has been done by iChiar & Tielensl ([2006) from the 
spectrum taken by the Infrared Space Observatory (ISO) to- 
wards WR 98a. This resulted in an extinction profile covering 
both the 10 and 20 /mi silicate fea ture. The spectral signature 
derived by IChiar & Tielensl d2006l) is much noisier than the 
10//m feature we used for our analysis above. Therefore, we 
chose to do the detailed analysis on the 10 //m feature and then 
check if the resulting dust composition and shape distribution 
is consistent with the full spectral profile. 

In order to check whether the composition and shape dis- 
tribution obtained from the fit to the 10 /mi feature is consistent 
with the observed 20 fim feature, we computed the expected 
extinction profile over the total spectral range for the dust com- 
position and shape distribution as derived from the 10 fim fea- 
ture. To account for possible errors in the intrinsic spectrum 
of the source used to extract the dust extinction profile from 
the observations, we added a smooth quadratic continuum to 
the computed extinction profile in order to fit the observed ex- 
tinction spectrum. We do not think that the exact choice of the 
continuum has a significant impact on our results. The results 
are shown in Fig. [8] It is clear from this figure that the spectra 
using (coated) porous GRF or DHS particles give a much bet- 
ter fit to the observed spectral profile than the spectrum using 
hom ogeneous spheres. Th is is consistent with the conclusion 
by IChiar & Tielensl d2006l) that perfect homogeneous spheres 
cannot account for the observed 20 fim feature. 

4. Implications 

4.1. Processing of dust in protoplanetary disks 

An important unsolved issue is that of the apparent difference 
in the composition of crystalline and amorphous silicates in 
protoplanetary disk material: silicate crystals are found to be 
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Spherical grains Porous GRF grains 


Coated GRF grains 


DHS (f m , x = 0.7) 


Crystalline SiC (Laor & Draine 1993) 


3.6% 


(3.7o-) 


4.2% (4.3cr) 


2.6% 


(5.3cr) 


Crystalline 0-SiC (Mutschke et al. 1999) 


1.0% 


(2.6cr) 


1.2% (2.4cr) 


0.3% 


(4.0cr) 


Amorphous SiC (Mutschke et al. 1999) 


5.8% 


(3.1cr) 


5.0% (2.9cr) 


4.1% 


(2.8cr) 



Table 3. The results of the analysis using various SiC measurements availab le in the literature. We conclude that the SiC in the 



ISM has a lattice structure closest resembled by that of lLaor & Drai ne ( 1993). 



1 — ' — ' — ' — 1 — I — r 
Spherical grains 




Figure 8. The extinction spectra of the best fit models com- 
puted over a longer wavelength r ange (black curves) and com- 
pared to the feature obtained bv lChiar & Tielensl (120061) (gray 
curves). 



magnesium rich and iron poor, while amorphous silicates have 
been anticipated to have more or less equal amounts of iron 
and magnesium. As it is usually assumed that the bulk of the 
crystalline silicates in disks are formed by thermal annealing 
of amorphous silicates, it is difficult to reconcile their differ- 
ent chemical composition. The composition of the crystalline 
component is a firm result since emission spectra from crys- 
talline silicates show very narrow emission features, the wave- 
length of which is a strong indicator of the magnesium content. 
Cometary and circumstellar crystalline silicates are almost pure 



magn esium silicates (see e.g. iJager et al.1 Il998t Fabian et al 
20011) . Such magnesium silicate cryst als are expected to form 
from direct gas phase condensation jGrossmanll 19721). U sing 
radial mixing and thermal processing of grains lGaill J2.004h pro- 
vides predictions of the composition and lattice structure of sil- 
icate grains in the region where in the solar system the comets 
have assembled. From these predictions it follows that the crys- 
talline silicates in this region should be only partly magnesium 
rich - the part which is formed from direct gas phase conden- 
sation very close to the star - and partly still reflect the com- 
position of the original amorphous silicates - the part which is 
formed by thermal annealing of amorphous silicates. The mag- 
nesium rich gas phase condensates are only expected to make 
up ~ 10% of the crystalline silicates in the regions where the 
comets have assembled and thus cannot explain the high mag- 
nesium fraction in cometary crystals. We have shown, however, 
that the interstellar amorphous silicates are already very rich 
in magnesium. Thus the amorphous material from which the 
crystalline silicates form is already magnesium rich. This pro- 
vides a very natural explanation for the composition of crys- 
talline silicates in com etary and protoplanetary dust. Indeed, 
Davoisne et al. I d2006l) show that when an amorphous silicate 
with Mg/(Fe + Mg) = 0.9 is annealed, the resulting crystalline 
silicates are almost purely magnesium rich. 



4.2. Depletion of elements along various lines of sight 

If the interstellar silicates are magnesium rich, the question 
arises: where is the iron? Unusual depletion of Si, Fe and Mg 
of lines of sight in t he Small Magellanic C loud have been re- 
ported dWeltv et all l200lk ISofia et all 120061) . In these lines of 
sight most Si and Mg are in the gas phase indicating that the sil- 
icate component of the interstellar dust is destroyed. However, 
the iron depletion is unchanged, i.e. most of the iron is still 
in the solid phase. This indicates that iron and magn esium are 



in different dust materials as suggested already by ISofia et al 



(2006). If the iron would be incorporated in the silicate lat- 
tice it would return to the gas phase along with the magne- 
sium and silicon when the silicate is destroyed. Our analysis 
indeed shows that the silicate is highly magnesium rich indi- 
cating that the iron must be in the form of refractory species 
such as metallic iron or iron oxide. In this scenario, when the 
silicate is destroyed, the iron can stay in the solid phase while 
the magnesium and silicon are returned to the gas phase. This 
thus provides a natural explanation of the unusual depletion 
patterns reported. 
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4.3. The nature of GEMS 

It has been suggested that GEMS (Glass with Embedded Metal 
and Sulphides), which are small subgrai ns found in many IDPs , 
are unprocessed interstellar dust grains ( Bradle y et al. , 1999b ). 
This raises the question: is the chemical composition of ISM 
grains that we have derived consistent with that of GEMS? 
The idea that GEMS are unprocessed interstellar grains is sup- 
ported by a similarity between the shape of the 10 fim feature 
of some GEM S and th at observed for the interstellar medium 
(Bradl ey et al, I fl999b). although also distinct differences can 
be identified dBradlev et all 11998). In addition, traces of expo- 
sure of GEMS to ionizing radi ation con s istent with an interstel- 
lar origin have been reported (IBradlevi Il994l) . Also, for some 
individual GE MS it is determined that they are of extr asolar 



origin (see e.g. lMessenger et all |2003e iFloss et al. , 2006). 



The amorphous (glassy) silicate in GEMS has a magnesium 
content of Mg/(Fe + Mg) « 0.9, consistent with our analysis of 
ISM grains. The iron found in GEMS is mainly in the form of 
inclusions of FeS and metallic Fe, only a minor fraction is in- 
side the silicate lattice. Among these components iron sulphide 
(FeS) is dominant. This is in contrast with studies of depletion 
of sulfur from the gas phase in the diffuse interstellar medium 
which show that the s ulfur in the ISM is predominantly in 
the gas phase (see e g. Ueda et all 120051 : ISavage & Sembachl 



19961 Snow & Witt , 1996h . Preliminary observations indicate 



that the iron sulphide inclusio ns in GEMS are lo cated prefer- 
ably at the edge of the grains ( Keller et all 20051) . This might 
indicate that metallic iron inclusions inside ISM grains reacted 
with sulfur to form the FeS inclusions in GEMS. If so, this pro- 
cess must have happened in the collapsing molecular cloud, or 
in the protoplanetary disk phase. 

The average GEMS sil icate has (Mg + Fe)/Si ~ 0.7 
dKeller & M essengel |2004 . much lower than we derive 
((Mg + Fe)/Si « 1.5, see Table HJ. This gives O/Si ~ 2.7 
in the GEMS silicates, implying they are, on average, almost 
pure pyroxene with a small amount of silica. This difference is 
consistent with the spectral position of the ma ximum absorp- 
tion re ported for GEMS which is around 9.3 yum (iBradley et all 
1999al) . while it is around 9.7 fj.m in the ISM. Considering the 
fact that decreasing the O/Si in the amorphous silicates tends 
to shift the maximum absorption towards shorter wavelengths 
(see Fig. HJ, we can understand this difference between the 
GEMS and ISM 10 /mi features in terms of a difference in aver- 
age stoichiometry. If we apply abundance constraints as deter- 
mined for GEMS in our fitting of the 10 //m feature we get no 
satisfying fit (x 2 > 90 for all particle shape models). The argu- 
ments above taken together suggest that GEMS are, in general, 
not unprocessed leftovers from the diffuse ISM. However, they 
may have formed in the collapsing molecular cloud from which 
the Solar system was formed. 



5. Conclusions 

We have shown that the interstellar silicate extinction spectrum 
can be fitted accurately using irregularly shaped particles. In 
previous studies the silicate extinction has been explained us- 
ing spherical amorphous silicate grains with equal amounts of 



iron and magnesium. However, interstellar silicates are most 
likely not perfect homogeneous spherical grains. Using irregu- 
larly shaped dust grains we find a reduced^ 2 that is almost six 
times smaller than that obtained using homogeneous spheres. 
Moreover, the spectrum of homogeneous spherical grains is in- 
consistent with the observed 20 yt/m extinction feature. Using 
irregular particle shapes, we analyze the lO/im extinction fea- 
ture in detail. The resulting dust composition and shape distri- 
bution is also found to be consistent with the observed 20 fim 
feature. 

To summarize we conclude the following concerning the 
composition of the interstellar silicates: 

- The amorphous silicates in the ISM are highly magnesium 
rich. We show that the average composition of the silicates 
corresponds to Mg/(Fe + Mg) > 0.9. 

- We confirm the findings of previous studies that the crys- 
tallinity of the silicates in the ISM is small, although we do 
find a small fraction of crystalline forsterite. We find that 
the crystallinity is approximately 1%. 

- We find a mass fraction of approximately 3% of silicon car- 
bide. The presence of SiC in the ISM has been expected 
since presolar SiC grains are frequently found in IDPs. We 
identify the presence of SiC grains in the ISM spectroscop- 
ically with more than 3<x confidence. 

- The average stoichiometry of the interstellar silicates is 
between olivine and pyroxene type silicates. We find that 
O/Si -3.5. 

Our analysis has several implications for the interpretation 
of silicate dust in various environments. These are: 

- The high magnesium content of amorphous silicates in the 
ISM provides a natural explanation for the high magne- 
sium content of crystalline silicates found in cometary and 
circumstellar grains. We argue that these magnesium rich 
crystals can be formed by simple thermal annealing of the 
magnesium rich amorphous silicates. 

- Unusual depletion patterns observed in several lines of 
sight through the interstellar medium show that the deple- 
tion of Fe is not correlated with that of Mg and Si. This is 
naturally explained when these components are not in the 
same dust materials, consistent with our findings. 

- We discuss the origin of GEMS in interplanetary dust parti- 
cles. These grains are suggested to be of interstellar origin. 
However, we show that the composition of the ISM sili- 
cates is not consistent with that of GEMS, suggesting that 
GEMS are not unprocessed remnants from the ISM. 

As a concluding remark we would like to point out that the 
use of homogeneous spherical particles to model silicate ex- 
tinction spectra gives results that deviate in important aspects 
from the values obtained using more realistic particle shapes. 
Since most dust grains in nature are not homogeneous spheres 
the use of homogeneous spherical particles should be avoided, 
even if they do provide a reasonable fit to the observations as 
it may yield spurious grain properties. In cases where mod- 
eling of realistically shaped particles is computationally too 
time consuming, the Distribution of Hollow Spheres provides 
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a good alternative. We have shown that this shape distribution 
simulates the properties of irregularly shaped particles, while 
computational effort is small. 
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Appendix A: Construction of porous Gaussian 
Random Field particles 

In order to construct the Gaussian Random Field particles 
(hereafter GRF particles), we employed the following algo- 
rithm. First we constructed a three-dimensional field with ran- 
dom numbers R^, (z, j, k = 1..M). From this we constructed a 
field djk according to 

G ijk =R ijk e-P dl (A.l) 
where p determines the size of the particles to be formed and 

d = Si - M/2) 2 + (j - M/2) 2 + (k - M/2) 2 (A.2) 

is the distance to the center of the field. The normalized 
Gaussian Random Field is obtained from Guk by taking its 
three-dimensional Fourier transform 

&ijk = T (G ijk ) I max T (G m ) . (A.3) 

v ' ijk v ' 

where T denotes the Fourier transform. From this normalized 
three-dimensional field we construct the particles by taking 
all points in space where > 0.5 to be inside the particle, 
while Qijk < 0.5 is taken to be outside the particle. In this way 
Qijk represents a field filled with a large number of particles. 
We have to separate these in order to get single particles. 

For the computation of the extinction spectra of these par- 
ticles we use the Discrete Dipole Approximation (DDA). In 
the DDA we approximate the particle by discretizing the parti- 
cle volume into small volume elements that each interact with 
the incoming light as dipoles. The immediately provides a 
grid for these dipoles. For the final DDA calculations we only 
selected those particles that contain between 2000 and 3500 
dipoles. 

In order to construct porous Gaussian Random Field parti- 
cles, we make vacuum inclusions in the particles by computing 
a second Gaussian Random Field, Q\;^ with a much smaller ef- 
fective size, p' . The points in space where this second field has 
values larger than a certain threshold are considered vacuum. 
By changing the threshold value and the size of the vacuum in- 
clusions, p', we can vary the degree of porosity and irregularity 
of the final particles. We have taken p' = p/50 and a threshold 
value of 0.1. 

The coated porous GRF particles are created in the same 
manner as the pure GRF particles but we take another threshold 
value below which the material is considered mantle material. 
In this way the outer regions of the grain become the mantle, 
while the inner regions of the grain (where the GRF has its 



highest values) consists of the core material. This threshold is 
chosen such that each particle is 'covered' with a mantle con- 
taining 30% of the material volume of the grain. 

We constructed an ensemble of 38 porous GRF particles 
and averaged the optical properties over this ensemble. Typical 
examples of the particles are shown in Fig.Q] An example of a 
coated porous GRF particle is shown in Fig. [2] 
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